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Abstract 

We present a next-to-leading order QCD calculation of the production rates 
of jets containing heavy quarks. This calculation is performed using the stan- 
dard Snowmass jet algorithm; it therefore allows a comparison with similar 
results known at next-to-leading order for generic jets. As an application, we 
present results for the inclusive transverse energy of charm and bottom jets 
at the Tevatron collider, with a complete study of the dependence on the jet 
cone-size and of the theoretical uncertainties. 
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1. Introduction 

The study of heavy quark production has provided some of the most interesting 
results in the physics of high energy hadronic collisions. The large available data sets 
of 6-hadrons have been used for precise measurements of spectroscopy and lifetimes, 
as well as for measurements of the production rates. The associated production of 
jets including 6-hadrons and W vector bosons has been used for the detection of the 
t-quark. Several signals for new physics, such as an intermediate mass Higgs or the 
supersymmetric partners of the top quark, could manifest themselves via the presence 
in the final state, among other things, of jets containing 6-quarks. A close study of the 
production properties of fe-jets in QCD is therefore an important phenomenological 
input for many of these searches. Calculations have been performed in the past for the 
production of 6-quarks at next-to- leading order (NLO) in QCD 0,0] • They have been 
used in comparisons with data measured at the SppS [|J and Tevatron ||] colliders. 
In this paper we present a calculation of the production rates of jets including heavy 
quarks, at NLO in QCD. The main difference between the study of a heavy quark and 
a heavy-quark jet is that in the former case one is interested in the momentum of the 
quark itself, regardless of the properties of the event in which the quark in embedded, 
while in the latter case one is interested in the properties of a jet containing one 
or more heavy quarks, regardless of the momentum fraction of the jet carried by 
the quark. A priori it is expected that variables such as the E T distribution of a 
heavy-quark jet should be described by a finite-order QCD calculation more precisely 
than the p T distribution of open quarks. This is because at high momentum large 
logarithms \og(p T /m) appear at any order in the perturbative expansion of the open 
quark p T distribution, due to the emission of hard collinear gluons. These logarithms 
need to be resummed, using techniques such as fragmentation functions f| . Collinear 
logarithms, on the other hand, are not present in the E T distribution of heavy-quark 
jets, since the jet E T does not depend on whether the energy is carried all by the quark 
or is shared among the quark and collinear gluons. The experimental measurement 
of the E T distribution of heavy-quark jets does not depend either on the knowledge of 
the heavy-quark fragmentation functions, contrary to the case of the p T distribution 
of open heavy quarks. Experimental systematics, such as the knowledge of decay 
branching ratios for heavy hadrons or of their decay spectra are also largely reduced. 

The contents of this paper are as follows: in section 2 we introduce our notation, 
definitions, and summarize the technique used in the calculation. Section 3 presents 
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our phenomenological results, with a complete discussion of the production rates and 
properties of charm and bottom jets at the Tevatron collider. Our conclusions can 
be found in section 4, while a thorough discussion of the technical details of the 
calculation is collected in the Appendix. 

2. Cross sections 

In this section, we briefly discuss the definition of the heavy-quark jet cross section 
in perturbative QCD. The interested reader will find a more detailed presentation in 
the Appendix. 

Thanks to the factorization theorem || , a cross section in hadronic collisions can 
be written in the following way 

da^ H2 \P 1 ,P 2 ) = ]T j dx 1 dx 2 f^\x 1 )f { b H2 \x 2 )da ab (x 1 P 1 ,x 2 P 2 ) , (2.1) 

ab 

where Hi and H 2 are the incoming hadrons, Pi and P 2 their momentum, and the sum 
runs over all the parton flavours which give a non-trivial contribution. The quantities 
d&ab are the subtracted partonic short-distance cross sections, which are calculable 
in perturbative QCD. Our aim is to evaluate these quantities for heavy-quark jet 
production. We will deal with da a b, which are directly related to Feynman diagrams; 
da a b are obtained from da ab by subtracting some suitable counterterms for the initial 
state collinear singularities ||. 

In perturbative QCD, a jet has to be defined in terms of unobservable partons. 
The definition can be rather freely chosen, the only constraint being that it must 
guarantee the infra-red safeness of the cross section. In the present paper we adopt 
the Snowmass convention J7JI, whereby particles are clustered in cones of radius R in 
the pseudorapidity-azimuthal angle plane. 

The calculation of the heavy-quark jet cross section is very similar to the one of 
the generic-jet cross section, but two important differences have to be stressed. By 
its very definition, a heavy-quark jet is not flavour-blind; we have to look for those 
jets containing a heavy flavour. Furthermore, the mass of the heavy flavour is acting 
as a cutoff against final state collinear divergences; therefore, we will not have to deal 
with many singular contributions, which are present in the generic-jet partonic cross 
section. 
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At the leading order, the heavy-quark jet cross section in hadronic collisions gets 
contributions from the following partonic processes 

9 + 9 -> Q + Q, (2.2) 
q + q -> Q + Q. (2.3) 

At the next-to-leading order, the cross section gets contributions from the radiative 
correction to the processes ( |2.2|) and (|2.3D , and from the tree amplitudes of the pro- 
cesses 

g + g Q + Q + g, (2.4) 
q + q -> Q + Q + 3, (2.5) 

<? + # -> Q + Q + g. (2.6) 

In the following, we will consider the case in which the heavy-quark jet contains the 
heavy quark. This is by no means restrictive; the heavy antiquark obviously can be 
treated in the same way. 

We write the leading-order contribution to the heavy-quark jet cross section as 

daS = S 2 M<$d<f> 2 diij, (2.7) 

where M^b * s the leading-order transition amplitude for the two-to-two process 
a + b — > Q + Q, squared and summed over spin and colour degrees of freedom, and 
divided by the flux and average factors; o?$ 2 is the two-body phase space for the QQ 
pair, 

diij = dE JT drijdn^- 2e) (2.8) 

is the measure over jet variables in 4 — 2e dimensions, and S 2 is the so-called mea- 
surement function, which defines the jet variables in terms of the partonic variables 
(although eq. (|2.7|) is completely general, the explicit form of the measurement func- 
tion depends upon the merging algorithm adopted: for a discussion on the use of the 
measurement function in jet physics, see refs.||§]). In eq. ( |2.8| ) we indicated with E JT 
the transverse energy of the jet, and with rjj its pseudorapidity; dQf ^ is the angular 
measure in 2 — 2e dimensions. The measurement function is actually a distribution, 
since the jet variables are implicitly defined as entries of some 5 distributions; this is 
the reason for the dfij on the RHS of eq. (|2.7|) . The explicit form of S 2 is given in 
the Appendix. 
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The next-to-leading order contribution is 

d<rS! = d<r%} + da, 



to 



(2.9) 



where 



da. 



ab 



S 2 M^d$2dtij 



(2.10) 



is the virtual part, due to the radiative corrections to the two-to-two processes, and 



is the real part, due to the two-to-three processes. Here <i$3 is the three-body phase 



account the fact that, in the final state, one additional parton is present (and therefore 
the jet definition has to be suitably modified). 

It is apparent that, at the leading order, the heavy-quark jet can only coincide 
with the heavy quark itself. Therefore, at the leading order, the heavy-quark jet cross 
section is identical to the open-heavy-quark cross section. At the next-to-leading or- 
der, the presence of a light parton in the final state enriches the kinematical structure. 
The heavy-quark jet can be the heavy quark, or it can contain the heavy quark and 
the light parton, or the heavy quark and the heavy antiquark. The heavy-quark 
jet cross section is therefore different from the open-heavy-quark one. Nevertheless, 
thanks to the non-zero mass of the quark, the structure of the singularities of the 
heavy-quark jet cross section is identical to the one of the open-heavy-quark cross 
section; a detailed proof of this statement is reported in the Appendix. We can then 
write the heavy-quark jet cross section at the next-to-leading order, 



da^=S,M^d^d^ 



(2.11) 



space for the QQ pair plus a light parton, and 5 3 is analogous to S 2 , but takes into 



da ab = da ( J + da ( J , 



(2.12) 




in eq. fl2.13| ). The Fortran code for da ° b p is available from the authors of ref. 0. 
For the present paper, we wrote a Fortran code for the evaluation of dA a b. 



3. The structure of heavy-quark jets at the Tevatron 



As an application of the formalism developed so far, we present in this section some 
results of interest for measurements at the Fermilab 1.8 TeV Tevatron pp Collider [TO]. 
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Figure 1: b-jet inclusive E T distribution in pp collisions at 1.8 TeV, for 
\tj\ < 1, R = 0.7 and fi F = /u R = /j,q (solid line). For comparison, we also 
show the open-quark inclusive E T distribution (dashed line). The component 
of the jet-like contribution due to jets containing both b and b is represented 
by the dotted line. 



We will consider jets containing either charm or bottom quarks. For these we will 
provide absolute predictions for the production rates as a function of the jet transverse 
energy and jet cone size R, and will explore the theoretical uncertainties associated 
to the choice of factorization (yU p ) and renormalization (/i R ) scales. Since we will be 
considering jets of energy much larger than the heavy-quark mass, the uncertainty 
associated to the mass values chosen for charm and bottom quarks is negligible, and 
will not be discussed. We will also study the fraction of heavy-quark jets relative to 
generic jets, and the fraction of 6-jets relative to c-jets. For this particular distribution, 
we will show that most of the uncertainties related to the choice of scales cancel out 
in the ratio, leaving a rather accurate NLO theoretical prediction. 

We consider jets produced within < 1, in order to simulate a realistic geo- 
metrical acceptance of the Tevatron detectors. We will use the parton distribution 



set MRSA' (TT|. Our default values of the parameters entering the calculations are: 
m c = 1.5 GeV, m b = 4.75 GeV, /j r 2 = /i F 2 = /i 2 = m 2 Q + E 2 T and R = 0.7. 
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Figure 2: c-jet inclusive E T distribution in pp collisions at 1.8 TeV, for 
\n\ < 1, R = 0.7 and fj, F = fj, R = fiQ (solid line). For comparison, we also 
show the open-quark inclusive E T distribution (dashed line). The component 
of the jet-like contribution due to jets containing both c and c is represented 
by the dotted line. 



We start from the absolute production rates. Figure [I] shows the prediction for 
the E T distribution of 6-jets at the Tevatron. For the purpose of illustration, and 
to provide a direct estimate of the effects calculated in this paper, we separate in 
the figure the contribution of the open-quark component. As indicated in the figure, 
the jet-like component, defined as the additional contribution to the open-quark one 
(the A term in eq. (|2.13| )), becomes dominant as soon as E T becomes larger than 
50 GeV. We also show the part of the jet-like component due to jets that include the 
bb pair (we will call these bb jets). The figure suggests that for this E T range and with 
R = 0.7 this is the dominant part of the jet-like component. This is consistent with 
the expectation that, for large enough E T and provided that the majority of the final- 
state generic jets are composed of primary gluons, heavy-quark jets are dominated 
by the process of gluon splitting, with the jet formed by the heavy-quark pair. As 
we will show later, the situation changes at higher E T values, where heavy quarks are 
mostly produced via the s-channel annihilation of light quarks. 

Figure || shows the same distributions for c-jets. Notice that the value of E T at 
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Figure 3: Left: relative contribution of the jet-like component in the b- 
jet inclusive E T distribution, for different values of the renormalization and 
factorization scales (/% = /x F = fi). Right: relative contribution of the jet- like 
component in the b-jet inclusive E T distribution, for various cone sizes R. 
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Figure 4: b-jet inclusive E T rate, as a function of the cone size R, for 
E T = 50 GeV (left) and E T = 100 GeV (right). 



which the jet-like component becomes dominant is smaller than in the case of 6-jets. 
Again this is in agreement with naive expectations. The relative probability of finding 
a heavy-quark pair inside a high-.E T gluon scales in fact like \og{E T / wig) [[TjJ. 

Figure ||] shows the fraction of the jet-like component of 6-jets versus E T , for various 
choices of factorization/renormalization scale and cone size. In the former case the 
dependence on scale variations is rather small, while the dependence upon the jet 



-8- 




Figure 5: Left: relative contribution of the jet-like component in the c- 
jet inclusive E T distribution, for different values of the renormalization and 
factorization scales (/x R = /x F = /j>). Right: relative contribution of the jet-like 
component in the c-jet inclusive E T distribution, for various cones sizes R. 
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Figure 6: c-jet inclusive E T rate, as a function of the cone size R, for 
E T = 50 GeV (left) and E T = 100 GeV (right). 



definition is more significant. Notice in particular that the E T value at which the 
jet-like component becomes dominant depends significantly on the cone size, being 
equal to 25, 50 and 100 GeV for R — 1, 0.7 and 0.4 respectively. The absolute 
rates at E T = 50 and 100 GeV as a function of the cone size are given in fig. As 
explained in the Appendix, the cross section at R = is well defined, and it is equal 
to the open-quark cross section. This should be contrasted with the case of generic 
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Figure 7: Scale dependence of the b-jet E T distribution (R = 0.4, solid lines) 
and of the open-quark inclusive E T distribution (dashed lines). 



jets, in which the cross section at R = is not well defined, being negative at any 
fixed order in perturbation theory ]13 . 



Similar results for c-jets are shown in figs. |5| and |6|. 

The strong scale dependence exhibited by the absolute rates is of the same size 
as the one present in the inclusive p T distribution of open heavy quarks. This scale 
dependence is usually attributed to the importance of the gluon splitting contribution. 
This process appears for the first time at 0(a^) and is therefore, strictly speaking, 
calculated with leading-order accuracy only. One expects therefore that in a regime 
in which the gluon splitting contribution is suppressed by the dynamics the scale 
dependence should be milder. We will show later that such a suppression takes place 
for high-energy heavy-quark jets. Figure [7] shows the scale dependence of the 6-jet 
cross section as a function of E T , for values up to 450 GeV. In the high-.E T region the 
scale dependence is indeed reduced to the value of 20% when the scale is varied in 
the range /io/2 < ji < 2/i , a result consistent with the limited scale dependence of 



the NLO inclusive-jet cross sections [[13 



In spite of the strong scale sensitivity at the smaller values of E T , it is reasonable 
to expect that the ratio of the b- and c-jet rates be a stable quantity under scale 



0.4 | — b to c ratio, R=0.4 
Solid: /2//i — 1 
Dashed: fj.//j. — 2 
Dotted: /j./fi Q =0.5 
0.2 | — Upper: b — quarks/c— quarks 
Lower: b— jets/c— jets 



0.0 




20 



40 60 
E T (GeV) 



80 



100 20 40 60 

E T (GeV) 



80 



100 



Figure 8: Ratio of the b-jet to c-jet E T distributions, for different values of 
the renormalization and factorization scales (n R = fi F = [i) for R = 0.4 (left) 
and R = 0.7 (right). For comparison we also show the ratio of open quark 
inclusive E T distributions of b- and c-quarks. 



variations. That this is indeed the case is shown in fig. |8|. 

Of direct interest for studies of heavy-flavour tagging and for searches of possible 
new physics is the fraction of heavy-quark jets relative to generic jets. This is also in 
principle the most straightforward measurement from the experimental point of view. 
We present in fig. © the ratio of the 6-jet and inclusive-jet E T distributions [O] (a 



similar plot for the c-jet fraction is shown in fig. [T(J). The inclusive-jet E T cross section 
used here was calculated with the JETRAD program |14|, using the same choices of 
parton densities and (/i R ,/i F ) as were adopted for the 6-jet and c-jet calculations. 
Contrary to the figures presented so far, which showed results for the heavy quark 
only (i.e. no antiquark contribution), we adopt for this one the prescription used in 
the definition of the data presented by CDF ||10||. The 6-jets are defined there as 



jets containing either a b or a b quark, jets containing both being counted only once. 
We will call these 6(6)-jets. This distribution can be obtained by subtracting the 
contribution of 66-jets from twice the total 6-jet rate. 

As shown in fig. 0, the normalization of the heavy-quark jet fraction depends 
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Figure 9: Ratio of the b(b)-jet to inclusive-jet E T distributions, for different 
choices of renormalization and factorization scales = fi F = fi) for R = 0.4 
(top) and R = 0.7 (bottom). The data points for R = 0.4 represent preliminary 
results from the CDF experiment [TdjJ , for which only the statistical uncertainty 
is shown. 



on the choice of scale as well as on the jet definition. In particular, the ambiguity 
induced by the change of scale is of 0(2). While this uncertainty prevents an accurate 
prediction of the heavy-quark jet fractions, it is important to point out that the choice 
of scale for this process is not independent of the scale chosen to predict the inclusive 
open-quark p T distributions. Since the data on both the bottom quark j|,[nj and the 
inclusive-jet E T spectra fIE|]rF|| strongly support the choice /z R , /i F ~ or smaller, 
we suggest that this is the best choice for the scales to be used in the prediction of 
the 6-jet fraction. 

It is interesting to notice that with this choice of scale there is good agreement 
between the theoretical prediction and the CDF data, at least in the case of 6-jets. 
Notice also that, while the data on inclusive 6-hadron distributions require the choice 
of even smaller scales (of the order of //o/4 ||15||), the measurement of the 6(6)-jet rates 
indicates that the choice \i = fio/2 is adequate. As for the large disagreement with 
the charm data, we have no significant comment to make. Hopefully, additional data 
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Figure 10: Ratio of the c(c)-jet to inclusive-jet E T distributions, for different 
choices of renormalization and factorization scales (fj, R =fj, F =fi) for R = 0.4 
(top) and R = 0.7 (bottom). The data points for R = 0.4 represent preliminary 
results from the CDF experiment [TdjJ , for which only the statistical uncertainty 
is shown. 



will soon be available, as well as estimates of the experimental systematics. Should 
the disagreement persist, this would indicate the presence of theoretical systematics 
not accounted for by the standard procedure of exploring the scale dependence of the 
rates. Notice that the largest contribution expected from higher-order perturbative 
corrections is given by the production of cc pairs from softer gluons emitted during 



the gluon shower evolution. However these effects have been estimated in ref. [[12 



and have been shown to be negligible at the energies of interest for the current mea- 
surements. 

To conclude this section, we discuss the behaviour of the b-jet production cross 
section at high E T . This item is interesting in view of the discrepancy reported 
by CDF [[RJ in the tail of the jet distribution. If this discrepancy could not be 
accommodated by new theoretical developments in QCD [18] or in the fitting of 
parton densities a study of the flavour composition of these high-energy jets 
could help in understanding the nature of the phenomenon. 
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Figure 11: Ratio of the b(b)-jet to inclusive-jet E T distributions for /% = 
/x F = /i/2 and with R = 0.4 (solid) or R = 0.7 (dashes). The data points are 
preliminary CDF data pffl , obtained with R = 0.4. 



Figure |ll| shows the 6(6)-jet fraction for E T values up to 450 GeV, for two different 
values of the cone size (R = 0.4 and R = 0.7) and at /i R = [i F = fio/2. Notice 
that while the fraction remains constant through most of the E T range, a rise is 
observed above 300 GeV. To better understand the origin of this rise, we present 



in fig. [T^a the separate contribution to the 6-jet cross section of the three possible 
initial states, gg, qq and qg. Notice that the qq contribution becomes dominant 
for E T > 250 GeV. Figures [l2|b-d show, for each individual channel, the separate 
contribution of the open-quark and bb-jet components. For E T large enough, the 
dominant component of the gg and qg channels is given by the 66-jet contribution, 
because of the gluon-splitting dominance. In the case of the qq channel, on the 
contrary, the 66-jet term is always suppressed, and most of the 6-jets are composed of 
a single b quark, often accompanied by a nearby gluon. We conclude that at high E T 
the dominant mechanism for the production of a 6-jet is the s-channel annihilation of 
light quarks. Since at high E T mass effects are negligible, 1/5 of the jets produced in s- 
channel annihilation are 6-jets. A simple LO calculation shows that the fraction of the 
two-jet rate due to s-channel light-quark annihilation is about 20% at E T = 450 GeV, 
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Figure 12: Initial state composition of the b-jet production processes, calcu- 
lated for fi R = fi F = n/2 and R = 0.4 (upper left). Different components of 
the production processes: gg — ► b-jet (upper right), qq — > b-jet (lower left) and 
qg — > (lower right). 



giving an overall 6-jet over inclusive-jet fraction of approximately 4%. This explains 
the rise of the 6-fraction at high E T , and provides a nice consistency check of our 
results. Notice also that while the probability that a gluon-jet will split into a bb 
pair grows at large E T faster that what predicted by the 0{o? s ) calculation Jl2|0, the 
fraction of primary gluons in the final state is so small that the overall effect on our 
predictions is negligible. 



This happens because of pairs emitted at higher-order from the gluons of the shower. 



-15- 



4. Conclusions 

We presented in this paper a calculation of the production of jets containing 
heavy quarks, at NLO in perturbative QCD. The techniques we employed represent 
a further elaboration of standard methods developed in recent years for the study of 
higher-order processes in hadronic collisions. As a phenomenological application of 
our results, we presented a detailed discussion of the rates and properties of charm and 
bottom jets produced at the Tevatron pp collider. We found that some distributions, 
such as the ratio of bottom to charm jets as a function of the jet transverse energy, 
are rather independent of the choice of renormalization and factorization scales. We 
presented results for the b- and c-jet fraction of inclusive jets, and discussed in detail 
the properties and composition of very high E T jets. We found that these vary sig- 
nificantly across the E T range measurable at the Tevatron. We provided predictions 
for the 6-jet fractions at high E T , to be compared with what could be measured by 
the Tevatron collider experiments. 
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APPENDIX A: Heavy-quark jet definition 

In this appendix we report the technical details that were not explicitly given in 
the previous sections. In the following, we will always indicate the kinematics of a 
partonic process as 

a( P i) + b(p 2 ) -> Q(fci) + g(fc 2 )+c(A;), (A.l) 
where c is the final state light parton present in the two-to-three processes. 
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We begin by considering the leading-order cross section, eq. fl2.7[) . As discussed 
in section 2, in this case the heavy-quark jet coincides with the heavy quark. The S 2 
measurement function formally states this obvious fact 

S 2 = 5{E JT - E 1t )8( Vj - rj^-^iipj - <p ± ). (A.2) 

Here we denote by E JT the transverse energy of the heavy-quark jet, and with t]j 
and ip j its pseudorapidity and azimuthal angle respectively. E 1T — k®/ cosh(r/ 1 ) is the 
transverse energy of the heavy quark; notice that the transverse energy is equal to 
the transverse momentum only in the case of massless particles. 

We then turn to the next-to- leading order contribution, eq. fl2.9Q . Due to eq. (|A.2| ), 
the heavy-quark jet cross section will be different from the open-heavy-quark one only 
in the real part, eq. ( 2.11 ). Using the merging algorithm of ref. |7j], we have 



= 8{E JT - E 1T )5{ Vj - Vi)^- 2 ^ - 

x 0^\u l -u\> g{E 1T ,E T ) S j9\\i0 1 - u 2 \ > g(E 1T ,E 2T ) S j 

+ 5{E JT - E 1T - E T )6 yr]j J 5 yipj J 

x 6 ^ \u 1 — oj\ < g(E 1T , E T )^ 
+ o[E JT — E 1T — E 2T )o j 6 \tpj J 

xe[\u t -u 2 \<g{E 1T ,E iT )\, (A.3) 



where 



and 



9(x,y)= X+ ( \ R (A.4) 
max(x, y) 



\ui -Uj\ = \l (r]i - r]j) 2 + (tpi - (fj) 2 . (A.5) 

In eq. ( |A.4|) , R is the usual jet-resolution parameter, which defines the cone size 
in the pseudorapidity-azimuthal angle plane. The quantities appearing in previous 
equations are expressed in the laboratory frame. 

It is well known that perturbatively calculated QCD cross sections have a diver- 
gent behaviour even after the ultraviolet renormalization, in those regions in which 
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a massless parton is soft or collinear to another massless parton. These singularities 
cancel after the virtual part, the real part and the collinear counterterms are added 
together. To perform this cancellation analytically is therefore important, in order 
to disentangle the various singular contributions to the cross section. To this end, 
the preliminary step in our formalism is to study the behaviour of the measurements 
functions in the soft and collinear regions. In the case at hand, the only non-trivial 
case is the one of the 1S3 function, since the two-to-three processes are the only ones 
in which a massless parton is emitted. It is straightforward to see that the following 
equations hold 

lim S 3 = S 2 , (A.6) 

k°—>-0 

lim (S3 = lim £3 = S2 ■ (A. 7) 

k — >ki k — >k2 

This properties guarantee the infra-red safeness of the heavy-quark jet cross section. 
Notice that, due to the fact that the heavy quark is massive, we do not have to care 
about the final state collinear emission, which is relevant, on the other hand, for the 
generic-jet cross section calculations. 

We now have to disentangle the singularities appearing in the heavy-quark jet 
cross section. For this purpose, we use the same method as was recently presented in 
ref. [0]. We parametrize the four- momentum of the outgoing light parton as 

k=^(h^-V 2 e T ,y), (A.8) 

where e T is a (2 — 2e)-dimensional unitary vector in the transverse momentum space, 
— 1 < y < 1 and < £ < £jr,; the collinear and soft limits are y— > ± 1 and 
respectively. With this parametrization, we have 



E T = -^/W . (A.9) 
Defining 

d& n \k) = " 1 , n , (A.10) 

v ; (27r) n - 1 2fc° ' v ; 



we get 



Vs 



-2c 



E-H^\ k ) = — * ( ^ ) r 1 -" (1 - y'Y^d^yd^ . (A.11 



2(2tt) 
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where dVlP' 2e ) is the angular measure in 2 — 2e dimensions. This form is suitable for 
disentangling the singular contributions to the real cross section. We use the identity 

(A.12) 



where 



f -2e 
>,cut 

2T 

{^Ucut, 



-1-e 



2e 
2 (?> 



*(0 (i - 2/ 2 ) 

— [5(1 -t/) +5(1 



T 



2e 



'log^ 



i + i/, 



(A.13) 
(A.14) 



The distributions in eqs. ( |A.13| ) and ( A.14 ) are defined as follows 

j) - 







/ 



-l l-j/ 

1 4 , /(y)-/(-iM-y-i + ycu f ) 
-i 1 + ?/ 



(A.15) 
(A.16) 
(A.17) 



where < £ cu t < 1 and < y cu t < 2 are arbitrary parameters, which can be freely 
chosen to improve the convergence of the results in the numerical computations. 
Writing the three-body phase space as 



rf$ 3 = E 2 T d^ 2 E~ 2 d^- 2t \k) 



(A.18) 



(by construction, d<& 2 is basically identical to the two-body phase space for the QQ 
pair introduced in eq. ( j2.7|) , except for the delta over the four-momentum, which in 
the present case reads 6(pt + p 2 — k\ — k 2 — k)), we can exploit eqs. ( |A.11|) and ( |A.12j ) 
to write eq. ( |2.11| ) as follows (we neglect 0(e) terms) 



da 



(r) 
ab 



S 3 M^E 2 T d$ 2 



+ S 3 M%E 2 d$ 2 



2(27r) 3 - 2e 
1 



'y/S s 



-2c 



-2c 



2(2 



7T 



|3-2e 



V(C,y)d^dydQ^d t x J 

K{£, y)dtldydn { - 2 - 2t) d^ J . (A. 19) 
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The measurement function S 3 in the first term on the RHS of eq. (|A.19|) can be 
substituted with S 2 . In fact, the 8(£) and 5(1 ± y) contained in the factor V(^,y) 
allow us to take the soft and collinear limits of that term, that is, to exploit eqs. ( |A.6|) 
and ( |A.7| ). We can also write 

S 3 = S 2 + S 3 (A.20) 

where 

S 3 = SiEjT-E^Sivj-vJS 1 - 2 *^-^) 

x - AMu^ -uj\ < g(E 1T ,E T )j - \u x - u 2 \ < g(E 1T ,E 2T ) 

+ 6(E JT - E 1T - E T )8 U - *\ + vEt ) 5^ U ^- + ^ 



E JT 7 V Ej_ 
x 0^\uj x -u\ < g(E 1Tl E T )^ 

+ o\E JT E 1T E 2T )o I Tjj )6 [<pj- 



Ejt 7 V E JT 

xe^\uj 1 -uj a \<g(E XT ,E 2T )^, (A.21) 

having used the identity 

9{a > b)6(c > d) = 1 - 9(a < b) - 6(c < d) + 9 (a < b)6(c < d) (A.22) 

in the first term on the RHS of eq. (|A.3| ). Notice that the first term on the RHS of 
eq. ( |A.21p should also have a contribution of the kind 

e( <g{E 1T ,E T )\d^-u 2 \ <g(E 1T ,E 2T )\ (A.23) 

Nevertheless, this quantity is always equal to zero, since otherwise it would be possible 
to merge in a jet all three particles in the final state. Equation (|A.19 ) therefore 
becomes 

-2e 



d"8 = S ^ E T d ^ 2 {2lf-^ \t) V&yHdydrt 2 - 2 ^ 
+ S 2 M$ E 2 T d$ 2 — \— (y^) K{& y)didyd^ 2 -^d^ 



2(2n] 



+ S 3 M^E 2 T d^ 2{2 ^ 3 _ 2e fffi K{tv)dZd V drt*-*>dvL 3 . (A.24) 
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Using again eqs. ( |A.12j ) and (|A.18|) we finally get 



da. 



(r) 
ab 



da 



(r,open) 
ab 



+ dA 



ab 



with 



da 



(r,open) 
ab 



S2M < £d$ 3 d(M 



(A.25) 



(A.26) 



dA 



ab 



S 3 M^E 2 T d^ 2 



-2c 



K(£, y)d^dydn {2 - 2e) dfij . (A.27) 



2(27r) 3 - 2e V 2 , 

Equation ( |A.26|) is just the real contribution to open-heavy-quark production. We 
have therefore 



da 



(open) 



da 



(0) 



'ab -^ab+daS+daZr™', (A.28) 
where the first two quantities on the RHS of this equation were given in eqs. (|2.7|) 
and Q2.10p , respectively. 

We observe that, although eq. (A. 20) has a general validity, in the case of generic- 
jet production it is totally useless, since eq. ( |A.27| ) would display divergences due to 
final state collinear emission. These divergences would be cancelled by corresponding 
ones in eq. ( A.26 ). On the other hand, in the heavy-quark jet case the mass of the 
heavy quark is acting as a cutoff against these divergences, and eq. ( |A.27| ) is finite and 
can be numerically integrated. In this sense, eq. ( A.27 ) is peculiar of the heavy-quark 
jet cross section. Notice finally that this term vanishes when R — > 0, that is, when no 
merging is performed; this is physically sensible, since in the absence of merging we 
expect the heavy-quark jet cross section to coincide with the open-heavy-quark one, 
as is formally stated by eq. ( |A.25| ). 

To integrate numerically eq. ( [A.27 ), we observe that in the collinear limits y — > ±1 
and in the soft limit £— >0 the quantity S3 vanishes. It follows that the subtractions at 
y = ±1 and £ = 0, implicitly contained in the factor 1Z(£,y), eq. ( |A.14j ), are actually 
immaterial, and we can therefore perform in eq. ( |A.27| ) the formal substitutions 

(A.29) 



(r, opera) 




l±y 
1 



Furthermore, since all divergences have been properly regulated, we can set e 
eq. (|A.27|) . After some algebra, we get 

S 3 M { 2 d^dEj^jdifj . 



dA 



ab 



(A.30) 
= in 

(A.31) 
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Using eqs. ( |A.28|) and ( |A.31| ) we get eq. (|2.13| ) and the explicit form for dA a t,. We 
finally observe that, since eq. ( |A.31| ) is free of divergences, it is suitable for numerical 
integration. 
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